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Abstract: We study a cosmological scenario in which inflation is preceded by a 
bounce. In this scenario, the primordial singularity, one of the major shortcom- 
ings of inflation, is replaced by a non-singular bounce, prior to which the universe 
undergoes a phase of contraction. Our starting point is the bouncing cosmology 
investigated in Falciano et al. (2008), which we complete by a detailed study of the 
transfer of cosmological perturbations through the bounce and a discussion of possi- 
ble observational effects of bouncing cosmologies. We focus on a symmetric bounce 
and compute the evolution of cosmological perturbations during the contracting, 
bouncing and inflationary phases. We derive an expression for the Mukhanov- Sasaki 
perturbation variable at the onset of the inflationary phase that follows the bounce. 
Rather than being in the Bunch-Davies vacuum, it is found to be in an excited state 
that depends on the time scale of the bounce. We then show that this induces os- 
cillations superimposed on the nearly scale-invariant primordial spectra for scalar 
and tensor perturbations. We discuss the effects of these oscillations in the cosmic 
microwave background and in the matter power spectrum. We propose a new way to 
indirectly measure the spatial curvature energy density parameter Q/c in the context 
of this model. 



Keywords: cosmology, bounce, inflation, CMB. 



Contents 



1. Introduction 2 

2. Background cosmology 4 

2.1 Preliminaries 4 

2.2 The de Sitter cosmological bounce 5 

2.3 Scalar field-driven cosmological bounce 6 

2.4 Scales and variables 8 

3. Scalar perturbations 9 

3.1 Scalar perturbation variables 9 

3.2 Time evolution of u 12 

3.3 Relating u and v 14 

3.4 Initial conditions for v 15 

3.5 Primordial scalar perturbation spectrum 16 

4. Tensor perturbations 19 

4.1 Tensor perturbation variable 19 

4.2 Primordial tensor perturbation spectrum 20 

4.3 Tensor-to-scalar ratio 21 

5. The characteristic time scale Ar] 22 

5.1 Inflationary fleld space dynamics 22 

5.2 Computing Arj 23 

6. CMB and matter power spectra 26 

6.1 Scalar anisotropy spectrum 27 

6.2 Matter density power spectrum 29 

7. Conclusions 30 

A. Time evolution of curvature perturbations 37 

A.l Numerical examples 37 

A. 2 Growth of perturbations in the contracting phase 39 



- 1 - 



1. Introduction 



The flatness, homogeneity, isotropy and monopole problems of standard Big Bang 
cosmology are easily resolved by invoking a primordial epoch of inflation [1, 2, 3, 4]. 
This framework also explains the origin of today's cosmic microwave background 
(CMB) anisotropics and large scale structure with compelling naturalness and ease. 
It traces these anisotropics back to the quantum fluctuations of the scalar field re- 
sponsible for inflation. As the universe grows exponentially, these perturbations are 
stretched outside the Hubble radius at early times, become classical, source the tem- 
perature anisotropics of the CMB and provide the seeds for structure formation. In 
addition, all pre-existing super-Hubble fluctuations are conveniently stretched out- 
side the Hubble radius today and can be safely ignored. The statistical properties of 
the CMB and of the large scale structure of the universe in the observable range are 
then directly related to the initial quantum state of sub-Hubble primordial fluctua- 
tions at the onset of inflation. The initial quantum state is unambiguously defined 
since the wavelength of all perturbations observable today was much shorter than 
the Hubble scale at the onset of inflation so that perturbations start their evolution 
in the adiabatic vacuum. With this initial state prescription, inflation yields a pri- 
mordial spectrum of density fluctuations, which, when evolved through the radiation 
and matter dominated epochs of standard cosmology, yields a CMB temperature 
anisotropy angular power spectrum and a matter power spectrum compatible with 
observations [5] , lending strong support to the inflationary hypothesis, and any model 
of the primordial universe lacking this unambiguous initial state prescription gener- 
ally suffers from a lack of predictability. 

Inflation however raises two new issues and leaves one unsolved. Firstly, the 
identity and origin of the scalar fleld that drives inflation remain, for lack of candi- 
dates within the Standard Model (SM) of particle physics, entirely unknown despite 
the fact that, in recent years, signiflcant progress has been made by implementing 
inflation into a high energy framework beyond the SM. For example, in string theory, 
although scalar degrees of freedom abound, the derivation of their potentials remains 
challenging and few candidate flelds are suitable for inflation without (a more or less 
signiflcant amount of) flnetuning. 

Secondly, the inflationary mechanism used to generate classical perturbations 
from the stretching and ampliflcation of quantum fluctuations leads to a conceptual 
issue of self-consistency. As spacetime expands during inflation, the physical wave- 
length of any given perturbation increases proportionally to the scale factor of the 
universe. Conversely, if we follow a perturbation backwards in time, its wavelength 
becomes shorter, possibly decreasing below the Planck length if the duration of 
inflation exceeds some 60 to 70 e-folds. However, close to the Planck scale, quantum 
gravity corrections become relevant. This is the well-known trans- Planckian problem 
of inflation. It has been studied both theoretically and numerically [6, 7, 8, 9, 10, 11] 
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by introducing modified dispersion relations or imposing non-vacuum initial con- 
ditions on the perturbations or by introducing a new cut-off scale of order of the 
Planck mass in the inflationary mechanism. Although they remain phenomenolog- 
ical, these models generically predict oscillations superimposed onto the standard 
(approximately scale invariant) spectrum of primordial fluctuations. In Refs. [7, 8] 
it was shown that these oscillations then carry through to the CMB temperature 
anisotropy and the matter density spectra. 

Thirdly, the usual models of inflation, even in the context of eternal infla- 
tion [12, 13], are not past complete, and inevitably contain a time-like singularity 
in the past at which the universe is of vanishing size. On the contrary, higher or- 
der curvature terms in the Einstein-Hilbert action [14, 15, 16, 17, 18, 19, 20, 21], 
non-minimal coupling of matter fields to gravity [22, 23, 24, 21] and the low energy 
effective actions of some string theories [25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35] do 
allow non- singular cosmological background solutions. In these, instead of shrink- 
ing to zero size, the universe experiences a bounce: starting out large, the universe 
undergoes a contracting phase until it reaches a finite minimal size, after which 
an expanding phase occurs. The existence of bouncing solutions in these high en- 
ergy effective theories warrants the study of simpler classical non-singular cosmo- 
logical models, in which background and perturbations are completely understood 
and tractable. This simpler class includes models with minimally coupled scalar 
fields [36, 37, 38, 39], scalar fields with non-standard kinetic terms [40] or in geome- 
tries that depart from the Friedmann-Lemaitre- Robertson- Walker (FLRW) geometry, 
i.e. Bianchi or Kantowski-Sachs spacetimes [41]. Another possibility is to consider 
unconventional forms of matter [42, 43, 44, 45], or a combination of radiation and 
scalar field matter [46]. For an extensive review, see e.g. Refs. [47, 48]. 

Although bouncing cosmologies have most often been discussed as alternatives 
to inflation, see e.g. Refs. [49, 50, 32, 51], here we shall see both the contracting and 
the non-singular bouncing phase as cosmological epochs that connect to inflation. In 
this paper, the cosmological singularity is replaced by a classical bounce which can 
entirely be described within General Relativity (GR), and after which inflation takes 
place as usual. Moreover, we will mostly focus, for computational simplicity, on a 
symmetric bounce, meaning that the rate of contraction (up to a minus sign) and 
its duration are the same as in the inflationary phase that follows. The transition to 
standard Big Bang cosmology after inflation through a phase of reheating remains 
unchanged from the standard case without bounce. 

To determine the viability of a non-singular cosmology, it is important to inves- 
tigate the stability and attractor properties of the background cosmology and the 
necessary level of fine-tuning required in order to obtain a bounce [36, 37, 52]. At the 
perturbative level, it is important to understand the transfer of cosmological pertur- 
bations through the bouncing phase, whether there can exist a (robust) prescription 
to choose the initial conditions for cosmological perturbations at some stage in their 
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evolution, and finally, whether such bouncing models have distinctive observational 
signatures. 

In this paper, we focus on the latter three issues, using the background cosmology 
of Ref . [36] . This model relies on a dynamical scalar field with a symmetric potential 
of the small field type and a positive spatial curvature term. This is the minimal 
possible setup in order to obtain a bounce but violate the least number of energy 
conditions [36, 53]. 

This paper is organized as follows. We first give a detailed overview of the 
background cosmology in Sec. 2, and an introductory discussion of the evolution of 
cosmological perturbations through the bouncing phase in Sec. 3.1. We then derive, 
in Sec. 3.2, the expression of the Bardeen variable at the onset of inflation at first 
order in slow roll by solving its mode equation prior to and after the bounce and 
by then applying the standard matching procedure at the bounce. We compute the 
corresponding form of the initial state of the mode functions of the Mukhanov- Sasaki 
variable describing scalar perturbations at the onset of inflation (Sees 3.3 and 3.4). 
This initial state is an excited state rather than the vacuum state, and it is shown that 
it can result in the presence of oscillations at either leading or subleading order in the 
primordial spectrum of scalar perturbations, see Sec. 3.5. Similar results are obtained 
for tensor perturbations in Sec. 4. While the amplitude of the oscillations depends 
only on the initial state of the cosmological perturbations prior to the bounce, the 
frequency of oscillation is given by a new time scale, which is discussed in Sec. 5. In 
Sec. 5, we also discuss a new possibility to indirectly measure the spatial curvature 
energy density parameter Qic in the context of the model proposed. Finally, in Sec. 6, 
we write down the form of the angular spectrum of CMB temperature anisotropics at 
low multipoles, commenting on the appearance of slow oscillatory features on large 
angular scales, and also provide numerical examples, using the CAMB code [54], of the 
angular power spectra of CMB temperature anisotropics and the spectrum of matter 
density fluctuations. In Appendix A, the evolution of the curvature perturbation in 
the contracting phase is discussed in some more detail, and some numerical examples 
of symmetric and asymmetric bouncing cosmologies are provided. 

2. Background cosmology 
2.1 Preliminaries 

We consider a homogeneous and isotropic FLRW universe with line element 



where a{ri) is the scale factor as a function of conformal time rj (where drj = dt/a, 
with t the cosmic time). The constant parameter /C describes the curvature of spatial 



ds^ = (?7) -dr]^ + 



dr^ 





1 -/Cr2 
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sections and can always be rescaled such that /C = 0, ±1. In standard inflation, one 
usually sets /C = 0. For the non-singular model studied in this paper, however, 
we shall require /C = +1 (see below). For a perfect fluid, with equation of state 
P = wp, where p and P and w are the energy density, pressure and equation of 
state parameter, respectively, the Friedmann equations derived from the Einstein 
equations and the conservation of the stress energy tensor read 

K{p + p) = ^{n^-n' + jc) , (2.2) 

k{p + 3P) = -^n' , (2.3) 
p' + 3-H (p + P) = , (2.4) 

where n = Sir/m'^^ with rrip^ the Planck mass, and where = a' /a is the (conformal) 
Hubble parameter, with a prime denoting a derivative with respect to t]. 

Consider a bounce occuring at conformal time t] = (or equivalently t = 0). 
At the bounce, the scale factor a reaches its (finite) minimum, Ob, while the Hubble 
parameter evaluated at the bounce, "Hb, is equal to zero and its derivative with 
respect to conformal time, T-L'y^, is positive. For a perfect fluid, the null and strong 
energy conditions derived from the Penrose-Hawking singularity theorems [55] read 

p + P>0, p + 3P>0. (2.5) 

It can be seen from Eq. (2.2) that the null energy condition is preserved at the 
bounce only if /C > "H^,. Recast in terms of cosmic time, this inequality reads 
IC/o? > i/b where Ht^ is the derivative of the (cosmic) Hubble parameter with respect 
to t evaluated at the bounce. Hereon, we set /C = +1 and this condition turns into a 
condition on the combination a^ifb [see Eq. (2.24)]. From Eq. (2.3), one finds that 
the strong energy condition must necessarily be violated in order to obtain a bounce. 
Finally, note that the spatial curvature term K,/o? dominates the dynamics at the 
time of the bounce. An inflationary phase is therefore required afterwards in order to 
drive the spatial curvature contribution to the energy content of the universe down 
to a level compatible with observations [56]. 

2.2 The de Sitter cosmological bounce 

Let us first consider the de Sitter solution. If w = —1, the energy content of the 
universe is in the form of a cosmological constant. For /C = 0, one has 

a(t)oce^*, (2.6) 

where the (cosmic) Hubble parameter H = d/a is a positive constant so that the 
scale factor of the universe increases monotonically. Re-expressing Eq. (2.6) using 
the conformal time coordinate, one has 

1 1 

a(r]) = -— and U = aH = — . (2.7) 
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For t — >■ oo, a(t) — )■ oo so that — )■ while if, formally, one lets the scale factor 
become arbitrarily small at early times, then rj — > — oo in the past. In this case, one 
thus has T] G ] — oo, [. On the other hand, in a /C = +1 de Sitter universe, one has 



a{t) = Ob cosh 



t 

Ob 



and 



tanh 



Ob 



(2i 



For t < 0, one then has H < and the universe contracts exponentially. At t = 0, 
the dimensionful scale factor a{t) reaches its minimum Ob, and the Hubble parameter 
vanishes. For t > 0, on the other hand, H > and the universe inflates. Re- written 
in conformal time, one finds 



a{ri) = Ob sec (r/) and ^^{v) = tan (77) , 
so that the range of r] for a de Sitter bounce is G ] — 7r/2, 7r/2 [. 



(2.9) 



2.3 Scalar field-driven cosmological bounce 

If instead of a simple cosmological constant with w = —1, the energy content of the 
universe is dominated by a scalar field 9? with energy density and pressure 



/2 



P 



2a2 



/2 



and 



2a? 



it follows from Eqs (2.2) and (2.3) that 



ft (p + P) = fi; 



12 



{n^ -H' + K) 



and 



k{p + 3P) = 
while from Eq. (2.4) one has 



-U' = -2KV{ip) 



12 



a?V{^)_ 



2H^' + a 



,d]/ 

dip 







(2.10) 



(2.11) 



(2.12) 



(2.13) 



As stated earlier, at the bounce, "Hb = and l-L'^ > and if /C > Ti'^, the null energy 
condition is verified. This also implies 93' G R so that ghost fields are avoided. Again, 
the strong energy condition is not preserved at the bounce and implies alV (v^b) > V^b 
in Eq. (2.12). If, as required in order to satisfy the observed flatness of spatial sections 
today, an inflationary phase lasting A^'i^f = 60 to 80 e-folds is to take place after the 
bouncing phase, one requires the usual slow roll conditions. 



aV {lp) > p'^, Up' > p" 



(2.14) 



to be verified for A^inf e-folds. It is therefore seen that the strong energy condition 
is consistently violated throughout both the bouncing and the inflationary phase. If 
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in addition, the cosmology is taken to be symmetric or close to symmetric, then a 
deflationary phase prior to the bounce is required. The strong energy condition is 
then similarly violated in part of the contracting phase as well. 

The cosmological evolution of the early universe is conveniently described in 
terms of a set of horizon flow functions defined by [57] 

H' d In I e,- 1 , , , , 

which during inflation are slowly varying functions of time with <^ 1 for all i. The 
first three horizon flow functions are related to the usual slow roll parameters e, 6 
and ^ by 

e = ei, 5 = ei--e2, ^=-6263, (2.16) 

so that the functions ei and 62 correspond to the inequalities in Eqs. (2.14) while the 
third is related to the time derivative of e and 6 through its definition in terms of ^. 
The field (p is in the slow-roll regime as long as ej <^ 1 for all i and inflation ends 
when ei ~ 1. 

In the remainder of the paper, we shall work at first order in slow-roll, meaning 
that the parameters themselves are considered constant and terms of O (e^) or 
higher are dropped. In the context of a symmetric bounce, the slow-roll contracting 
(deflationary) phase is obtained by imposing the same conditions on the ej's but with 

n<o. 

In the standard inflationary picture, the horizon flow functions express the devi- 
ation from a pure de Sitter expansion in regions where the spatial curvature term is 
negligible. Here in fact, both the deflationary and inflationary background evolutions 
away from the bouncing phase itself can be fully described in terms of the e^'s. In the 
bouncing phase itself, where spatial curvature is important, the horizon flow func- 
tions are ill-defined^, but it is also possible to write down a generalized form of the 
scale factor describing deviations from a purely de Sitter bounce [38]. To describe 
the scale factor of a quasi de Sitter bouncing universe, we simply rewrite Eq. (2.9) 
in the generalized form 

a = ab cosh (ut) , (2-17) 

where u is some dimensionful scale characterizing the deviation from a de Sitter 
bounce for which u = l/a^,. In conformal time, one has 

a (77) = Ob sec ( — ) , (2.18) 

\VcJ 



^In fact, regardless of spatial curvature, in a bouncing cosmology with a bounce occuring at 
say t = 0, the horizon flow functions are well-defined only in the separate contracting (t < 0) and 
expanding (t > 0) domains but not globally since 'H{N) is not bijective throughout and since 
diverges at the bounce. 
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where ?7c = 1/ (wob) is a dimensionless conformal time scale. For a symmetric bounce 
(where contracting and expanding phases last for the same number of e-foldings A^) 
a Taylor expansion of Eqs. (2.11-2.13) and (2.18) around r] = yields [38, 36] 



2/C 



,2 ' 



(2.19) 



while the first Friedmann equation gives 

,9 QIC 

<^b = — 



- 2alV{(pi 



(2.20) 



where everywhere in these expressions the subscript "b" denotes quantities evaluated 
at the bounce. Eqs. (2.19) and (2.20) imply that 



UJ 



2al 



(2.21) 



If <y9(j = 0, a; = l/flb and the de Sitter solution is recovered. The null energy condition 
(2.5) is preserved provided ip'^ > 0, so that u < a^^. 

Let us discard the possibility that a bounce occurs when ip[^ = 0, in which case, 
for any bounce, the potential for the scalar perturbations diverges at the bounce, 
see Eq. (3.4) and Eq. (3.7). Then, once again under the assumption of a symmetric 
bounce, a bouncing solution is obtained only if [38, 36] 



(fb = 



dV 

dip 



0, 



dV 



dip^ 



< 0. 



(2.22) 



Note that these conditions can be expected to hold approximately for a (slightly) 
asymmetric bounce as well. The simplest renormalizable potential that satisfies the 
conditions (2.22) and is bounded from below reads 



V {<p) = 



(2.23) 



with M and fi a priori free parameters. Given this form of potential, the phase of 
inflation that follows the bounce falls into the class of small field models, for which 
the accelerated expansion sets in near the origin (p = 0. Requiring that both (p'^ and 
u be real further yields the constraint 



{T/k)M- 



2 < T < 3. 



(2.24) 



2.4 Scales and variables 



Having established the background cosmology, let us now briefly recall the relevant 
observable range of scales and orders of magnitude of the parameters of the potential 
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that will be used in the numerical work and observational predictions in Sees 5 and 
6, and in Appendix A. The Hubble scale and curvature radius today are given by 

^H = ^ and 4 = ^. (2.25) 

Given that the Hubble scale today is Hq = lOO/ikms"^ Mpc~^ with the reduced 
Hubble parameter h ^ 0.7, one has in = 3/i~^ Gpc. Furthermore, for /C = +1, given 
that \Q/c\ < 10~^ [5] and using = /C/(ao-f^o)^! one finds 

-f < 10"% (2.26) 

so that ao > 30h~^ Gpc. 

Although part of the motivation for studying non-singular cosmologies is the 
presence of bounce-inducing corrections to GR (or other bounce-inducing mecha- 
nisms) arising in some high energy theories, in this work, we shall remain in the do- 
main of applicability of GR and quantum field theory. We therefore set ~ 10^£p\. 
From the first Friedmann equation, one can see that this corresponds to an order of 
magnitude M ~ O{10~^mp^) for the mass parameter of V{ip), consistent with small 
field infiationary model values for M [58], see Eq. (2.28) but also Eq. (2.24) where the 
range of allowed M parameter values is given. Using these values for ao and ab, the 
number of e-foldings of expansion since the bounce then is = ln(ao/ab) > 125. The 
entire range of observable scales today, given in terms of the corresponding physical 
wavenumber fcphys, is 

10-3 hUpc < A;phys < 10=^ hUpc. (2.27) 

For ao = SO/i^^Gpc, the range of comoving wavenumbers is then 10^ ^ k < 10^. 
Finally, in order to obtain 60 to 80 e-foldings of infiation, the parameters of the 
potential (2.23) must be of the order of the following small field infiation fiducial 
values [58] 

M ~ C (lO-^mpJ , /i ~ C (mpj . (2.28) 



3. Scalar perturbations 

3.1 Scalar perturbation variables 

We shall now study the evolution of cosmological perturbations around the homo- 
geneous and isotropic background cosmology discussed in the previous section. We 
mainly consider the scalar part of the metric perturbations (given by the gauge invari- 
ant gravitational Bardeen potential $) in the presence of scalar field perturbations 
S(p, but also give, although in less detail, the result obtained for tensor perturbations. 
In longitudinal gauge, the scalar part of the perturbed metric reads [59] 

ds^ = (r/) [- (1 + 2$) dr/2 + (1 - 2$) -fffdx'dx^] , (3.1) 
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where '^l^ is the background metric of the spatial sections. Expressing the scalar 
parts of the perturbed Einstein equations in terms of $ and the density perturba- 
tion variable S^p and combining them in the appropriate way yields a second order 
differential equation for the modes $fc [59] 



II 



= 0. (3.2) 



and 



The wavenumber k is the eigenvalue of the Laplace-Beltrami operator on positively 
curved spatial sections. It is therefore a function of an integer n, and is given by 
k = ^yn{n + 2) with n = and n = 1 corresponding to gauge modes [60]. As is 
well-known, it is possible to define the variable u related to $ through the relation [59] 

$ = |(p + P)V2. = ^^«, (3.3) 
in terms of which the mode equation (3.2) reduces to 

< + [e - Vuiv)] = 0, (3.4) 

where 

Km4h-3K(1-c9. (I) (3.5) 

The quantity Cg can in some regimes be interpreted as the velocity of sound. The 
explicit form of Vu{r]) reads 

( //\ 2 /// 
^) -^-n' + AJC, (3.7) 
if' J if' 

which is manifestly regular when "H = 0. Another possible choice of variable is given 
by the generalized form of the Mukhanov-Sasaki variable, valid for arbitrary /C and 
which we shall denote v. The modes of the variable v are related to in Newtonian 
gauge, by the relation [38, 61, 62] 

where ^ 

X^ = l-3/C^^. (3.9) 

The gauge-invariant combination of the modes <I>fc and 5(pk defining the generalized 
Mukhanov-Sasaki variable in Eq. (3.8) yields an equation of motion for which is 
reminiscent of that for Uk , 

v'^+[e -V,ir])]vk = 0, (3.10) 
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where 



V.iri) = ^ + 3}C{1 - cl) and 4 = a-^. (3.11) 

The variable v has the additional property that it appears naturally as the canonical 
variable in the total action for cosmological perturbations from which Eq. (3.10) is 
derived and, as such, should be the starting point for quantization and for the choice 
of the quantum initial conditions seeding the growth of cosmological structures. It 
can however be seen from Eqs. (3.8-3.10) that using the combination defining v in 
Eq. (3.8) to construct the relevant variable for the study of the evolution of cosmo- 
logical perturbations is inappropriate in the present context. The use of v across 
the bounce is forbidden because the variable itself, Cg and the term z'l/zk are either 
ill-defined or divergent when 71 = [38, 63, 64, 45]. Furthermore, the potential Vy{ri) 
is fc-dependent so that defining an asymptotically flat spacetime for high frequency 
modes is ambiguous. From the expression for Xk, and setting /C = +1, there would 
also appear to be divergences whenever Cg = 1 — k^/3 [38]. This equality can be 
expressed as 

^2 in" 

where S is the second slow roll parameter defined in terms of the horizon flow func- 
tions in Eq. (2.16). In the context of small field models and away from 71 = 0, 
the slow roll parameter 6 verifies 6 < and \6\ <^ 1, so that the condition (3.12) 
is clearly not a problem. It is therefore possible to retain the generalized forms Vk 
and Zk far from "H = 0. In such regions, the modes Vk are everywhere well-defined. 
Furthermore, in the regions where both /C/a^ <^ 1 and where slow-roll applies, the 
cosmic Hubble parameter H is slowly varying and different from zero. One may 
therefore neglect the terms proportional to /C, and use the approximation Xk — ^■ 
The expressions for the modes Vk and for the /c-dependent function Zk then reduce 
to their standard (/C = 0) forms Vk and z, with z now /c-independent. A quantum 
interpretation of v is then possible and the standard quantization procedure on Vk 
can in principle be achieved in the usual way [59]. Furthermore, in these regions, the 
modes of the variables u and v may be related to one another in the standard way. 

t;. = (^y^^ (y)' and k\k = -z['^)'. (3.13) 

Using these relations, it is therefore possible to prepare an initial quantum state by 
choosing the mode functions Vk of the variable v in the regions 77 < where modes 
are sub-Hubble and spatial curvature is negligible but study the evolution of cosmo- 
logical perturbations through the entire bouncing phase using u. The resulting form 
of u after the bounce can then be used to determine the form of the modes Vk at the 
onset of inflation and then compute the primordial spectrum of scalar perturbations. 
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In the following sections, we study, analytically, the evolution of the perturba- 
tions using u after setting the initial state of perturbations using v. In Appendix A, 
the time evolution of perturbations and the issue of the growth of perturbations in 
contracting spacetimes is also addressed [65, 66, 67]. 

3.2 Time evolution of u 

Hereon, we restrict the analysis to the case of a symmetric bounce and in this section, 
we consider the evolution of the variable u in the contracting phase, at the bounce 
and in the inflationary phase. As shown in Ref. [36], the potential (3.7) is very well 
approximated by neglecting the curvature term /C/'H^, neglecting the time depen- 
dence of the horizon flow functions, and joining a slow-roll exponentially contracting 
phase directly to a slow-roll inflationary phase. This is true because in the range 
of wavenumbers of observational interest, the spatial curvature term and the small 
amplitude features in Vu{rj) near 77 = are always negligible compared to k'^. Note 
that, as discussed in the previous section, this approximation does not hold true in 
the case of the variable v. Let us also point out that while it can be neglected when 
studying the perturbation u, the contribution of /C, at the level of the background 
cosmology, is crucial to obtain the bounce. 

We now proceed as in Ref. [36]. We neglect the spatial curvature term, in which 
case, at first order in slow roll, Eq. (3.7) simplifies to 

Vuiv) = [e, + I) = ^ (ei + I) ' (3-14) 
where we have used the expression 

n^v) = = (3.15) 

r]-r]± x± 

In the last step, we defined x± = \f] — ^^±1, where the "— " or the sign is chosen 
according to whether one considers times 77 before or after the bounce at 77 = 0. Up 
to first order in slow roll, the equation of motion for u therefore simply reads 



u" + 



u = 0, (3.16) 



where the subscript k on the mode function u is implicit. Note that Eq. (3.16) applies 
for symmetric bounces only because for an asymmetric bounce, the values of ei and 
€2 cannot be assumed to be equal before and after the bounce. In Appendix A.l, we 
give numerical examples showing that in cosmologies with strong departures from 
symmetry, additional features appear in the potential for u. 

We now briefly discuss the meaning of the conformal time parameters ?7_ and 
?7+. The expressions for 'H± used in (3.15) are the analogs of the familiar expression 
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used in standard inflation where one has, to first order in slow roll, 

niv) = -—, (3.17) 
V 

with rj — oo, 0[. Here, instead, the conformal time intervals of the contracting and 
expanding phases in Eq. (3.15) are ]77-,0] and [0,?7+[ respectively, so that rj^ and 
?7+ are the asymptotic values of the conformal time for early and late cosmic times, 
which we shall label t_ and respectively. For instance, in the de Sitter bounce 
discussed in Sec. 2.2, one has 'r]± = ±7r/2. 

We now consider Eq. (3.16) separately in the two eras before and after the 
bounce, where the solutions for u take the form of Hankel functions of the first and 
second kind hI^^ and Hi^\ We can write the solutions explicitly as 

u-{r]) = [U^{k)Hl^\kx^) + U^{k)Hl^\kx^)] , (3.18) 

u+{r]) = [Utik)Hl^\kx+) + U^ik)Hl^\kx+)] , (3.19) 

where the order u of the Hankel functions (equal before and after the bounce) is 
given by 

^=^ + ei + |, (3.20) 

and where as before 

x_ = |?7 — ?7_| = ?7 — ?7_ , x+ = |?7 — ?7+| = ?7+ — ?7 . (3.21) 

It is important to realize that the arguments of the Hankel functions are different so 
that the expressions given in Eqs. (3.18) and (3.19) are distinct. These two solutions 
(and their derivatives) can be matched at the time of the bounce 77 = using the 
standard procedure. Doing so and expanding the Hankel functions for large values 
of \kx±\ (that is, for 77 ~ 0), one obtains 

U+ = {ak + i) e-^('=^^-"'^), (3.22) 
U+ = ((Tfc - i) e^^'^^^— '^), (3.23) 

where we retain only the leading order terms and where we have defined the param- 
eter 

and rewritten the conformal time difference as 

Ar] = ?7+ — ?7_ . (3.25) 

Note that the expression for ak is first order in slow-roll. This quantity will be seen 
in Sec. 3.5 to correspond to a slight tilt and amplitude modification of the standard 
form of primordial spectrum of scalar fluctuations. The quantity Ar] is the conformal 
cosmological bouncing scale which we shall discuss in some more detail in Sec. 5. We 
shall see that it is interpretable as the angular frequency of oscillations induced by the 
bouncing phase in the primordial power spectrum of scalar and tensor perturbations. 
In the case of a pure de Sitter bounce, Arj = ir. 
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3.3 Relating u and v 

In this section, we compute the form of the Mukhanov-Sasaki variable at the onset 
of inflation in terms of its form prior to the bounce. Once determined, the expression 
for the Mukhanov-Sasaki variable in the expanding phase can be used in the standard 
way to compute the primordial spectrum of scalar perturbations. 

Let us recall that when expressed in terms of m, the time evolution of scalar 
perturbations is continuous across the bounce. The coefficients of the two linearly 
independent solutions for can then be expressed in terms of those for . This 
calculation was the subject of Sec. 3.2. On the other hand, the time evolution in 
terms of the Mukhanov-Sasaki variable v is singular at the bounce and in fact v itself 
is ill-defined at the bounce (see Sec. 3.1). As a result, a calculation similar to the one 
performed in Sec. 3.2 is not possible in the case of v. One may nevertheless easily 
determine the functional form of the Mukhanov-Sasaki variable in the defiationary 
and infiationary phases on either side of the bouncing phase because these are regions 
where /C/a^ -C 1 and where slow roll applies so that -u ^ w and Vr,(?7) ~ V^irf). In 
these regions, the equation of motion for the mode functions v is then given, at first 
order in slow roll, by the standard expression 



v" + 



1 / 362 

' 2 + 3ei + 



xl V ' ' 2 



0, (3.26) 



where the subscript on f is implicit. The solutions of Eq. (3.26) are given by 

v~{r]) = ^/k^ [V^~Hl^\kx^) + V-fH^^\kx.)] , (3.27) 
v+{r]) = [V+H^'\kx+) + V+Hl'\kx+)] , (3.28) 

where g = u + 1 with u defined in Eq. (3.20). From the relations in Eqs. (3.13) 
between u and v, and using Eqs. (3.18), (3.19), (3.27) and (3.28), one finds 

^^ = ±^?)''^"V,^ (3.29) 



.3 

so that the coefficients of v~ and are simply related by 

Vi^ik) = Viik) (ak + t) e-*^'^^^--^) , (3.30) 
V+{k) = Vf{k) ((Tfc - t) e^C^^^^-f) . (3.31) 

We stress once more that although analogous to Eqs. (3.22) and (3.23), the expres- 
sions (3.30) and (3.31) cannot be obtained using the procedure followed to relate 
and u~ because, as discussed above, the potential Vy^r]) in Eq. (3.11) is singular at 
the bounce, and because v is ill-defined at the bounce. 

We shall now require that the canonical variable v be quantized. In regions where 
modes are sub-Hubble and spatial curvature is negligible, the standard quantization 
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procedure of the canonical variable v can be applied. The consistency of the temporal 
mode functions v with the commutation relations required for quantization is ensured 
provided v satisfies the Wronskian condition [59] 

(v)' V* - {v*y V = t . (3.32) 

This imposes 

|V;±|2_|y±|2 = ^^A;-i (3.33) 
on the coefficients of the Hankel functions in (3.27) and (3.28). 
3.4 Initial conditions for v 

Unlike for inflation, there exists no unique prescription to fix the initial conditions 
for the cosmological perturbations at the onset of the defiationary phase in a bounc- 
ing cosmology. This is because, rather than starting sub-Hubble in their vacuum 
state, modes that enter the contracting phase have a prior history. For instance, and 
although this is by no means the only possibility, they may have started sub-Hubble 
during a radiation- or matter-dominated contracting era connected to the deflation- 
ary phase that preceeds the bounce. At the onset of inflation, therefore, modes 
cannot be expected to be in their vacuum state. In this paper, we shall nevertheless 
assume, for deflniteness, that V{' and can be parameterized as 

= ^ ?i k-'^^'e'"' and V," = ^^2 k-^'^e^'^ , (3.34) 

where ^1, ^2 ^ IR and positive, a and (5 are numbers while 9i and 62 are phase angles. 
Given Eq. (3.34), one flnds from Eq. (3.33) that 

c^lk-^ = c^lk-'' - k~^ . (3.35) 

In order to satisfy Eq. (3.35), one has to choose a = P = 1 and hence ^2 — ^1 ^ ^■ 
Dropping the subscript on so that ^1 — )■ one flnally has 

V- = ^ 1^1 k-'/'e''^ and = ^ \l - k-^'^e^'^ , (3.36) 

where we have absorbed the choice of sign in going from | Vi \ ^ and | | ^ to and 
V2 into the phases 61 and 62- Additionally, it can be checked using Eqs. (3.30) and 
(3.31) that for 77 > the normalization condition is verifled up to flrst order in slow 
roll, 

\V,^?-W? = -\k^'[l + aD . (3.37) 

where cx^ is second order in slow roll, as can be checked from (3.24). 

The initial conditions for v~ at times 77 < have a simple interpretation. They 
represent, through the phases 61 and 62 and through the amplitude parameter a 
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deviation from the Bunch-Davies vacuum, which can be recovered provided ^ = 1 
and 6*1 = (7r/2)f) + 7r/4. Furthermore, provided one does choose the adiabatic vacuum 
prescription, it can be seen from Eqs (3.30) and (3.31) that is given by a single 
mode function. Consequently, in this case and at first order in slow roll, we can 
expect to recover the standard form of the primordial spectrum of scalar fluctuations 
(see Sec. 3.5). 

A comment on the order of magnitude of ^ is also in order. It is well known that 
cosmological perturbations on super-Hubble scales generically grow in contracting 
spacetimes [65, 66, 67]. This is because while in an expanding phase, a perturbation 
can be split into a constant and a decaying mode, in a contracting phase, a pertur- 
bation is split into a constant and a growing mode. The amplitude of a perturbation 
with wavenumber k, which we parameterize by ^, is therefore dependent on the rate 
and the duration of the contracting phase. As discussed in Appendix A. 2, the am- 
plitude of perturbations can be made to deviate only slightly from the amplitude of 
vacuum fluctuations provided the contracting phase is either slow, short, or both slow 
and short. In the case at hand, the contracting phase is a fast and short deflationary 
phase (~ 60 e-folds in a short conformal time interval) during which super-Hubble 
modes grow only very little. Furthermore, fluctuations can be assumed to remain 
small in any contracting phase preceding deflation, provided the rate of contraction 
is small. Finally, since <^ turns up in the expression for the amplitude of the pri- 
mordial, CMB and matter power spectra, that it remains close to 1 is evidently an 
observational necessity for the viability of the model. We shall therefore assume that 
^ ~ 1 when we consider cosmological observables in Sec. 6. 



3.5 Primordial scalar perturbation spectrum 

We shall now use the late time (super-Hubble) behaviour of v'^ to compute the 
primordial spectrum of scalar perturbations in the expanding phase, P^, where, in 
the subscript, the variable ( = v/z is the curvature perturbation in the comoving 
gauge discussed in the Appendix. The spectrum is given by the usual expression 



2^ 



2^ 



(3.38) 



and on super-Hubble scales, is given in terms of the asymptotic forms of the 
Hankel functions for kx+ <^ 1, 



Vk^\v,^{k) 



+V,+ {k) 



kx^ 



T{1 + q) V 2 



r(l — q) sin TTQ \ 2 



kxA 



r(l — q) sinvTf) \ 2 

' kx ^ ^ 



(3.39) 
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Given that, as long as slow roll is not violated, g > 0, we may keep only the second 
term in each line of (3.39). From the growing modes of (3.39) and from 



where, since /C/a^ ^ 1 we have taken x — ^ind where we have used the definition 
of the first horizon flow function, we find 

n = 2 (kx^f^'^ |V+ - y^+p. (3.41) 

m^^ TT a'ei [r(l - g) 

It is convenient to split the explicit expression for V( into a first standard (tilt) part 
P^**^ and a second part V^^'^, which will be seen to be oscillatory in k. Let us thus 
write 

= Vf^ X . (3.42) 
The standard part of the spectrum is given by 

^ 2ml^ahi[T{l-g)sm{ng)f^ 
which, up to first order in slow roll, reads 

pstd ^ 1^ _ ^ ^ (j^^ ^ ^2ei + 62) Hkx+)]} , (3.43) 

where C = 7e + ln2 — 2 and 7e is Euler's constant. Eq. (3.43) is the standard 
inflationary result for P^. The oscillatory part of the spectrum reads 

Vr = -k\V^^-V2^\\ (3.44) 
From (3.30) and (3.31), one has 

V+ -V^ = A {\V^\e'^''-^^ - \V^\e'^''+'^^) , (3.45) 

where 

1 /2 

A=[l + a^l) , = [(2ei + 62)0-^ ^ - TT^] - V9, v9 = arctana^^ (3.46) 

with (Tfc defined by Eq.(3.24) and Ar/ by Eq. (3.25). The expression in Eq. (3.45) can 
also be rewritten as 

V+ -V^ = A[\V^\'^ + \V^\^ - 2|yf ||V^2"I cos (^1 -62- 20)] ^/^e'*, (3.47) 
where 

ril/f |sin(ei -0) - |K^|sin(e2 + 0)l ,0 .on 

m = arctan ^-^ . (3.48) 



\V{- 


sin {61 — (f)) — 


V2^ 


sin {62 + 0) 




cos (^1 — 0) — 


V2~ 


cos [62 + 0) 
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Using these expressions, the oscillatory contribution to the spectrum V^^'^ reads 



pose _ 1 ^2 ^ r|y-|2 ^ m-|2 _ 2 1 1 1 1 COS {01 -02- 20)1 . (3.49) 



Up to this point, our calculations were exact to first order in slow-roll. We now make 
an approximation to simplify the cosine term in Eq. (3.49). We first rewrite it as 

cos (6*1 — 02 — 2(p) = cos (6*1 — 02 — 2k At] + 2nu) cos (2 arctancr^"'^) + 

sin (6^1 -02- 2k Ar] + 2'itu) sin (2 arctan a^^) . (3.50) 

Because ^ 1, and as k grows, the cosine and sine terms that have the arctangent 
as their argument rapidly go to —1 and respectively, so that one is simply left with 
just the term cos {0i — 02 — 2k At] + 27ri/) in the first line of Eq. (3.50). As for A^, 
one can easily see from Eq. (3.46) that it is equal to one at first order in slow-roll. 
Furthermore, expanding V^^'^ at that order and using Eq. (3.36) one obtains 

pose ^ ^2 ^ 1^ _ ^2| _ 2^|1 _ ^2|l/2 [cos(2fcAr7 + ^1 - 02) 

-7r(2ei + ea) sm{2kAr] + 0^ - ^2)] • (3.51) 

The V^^'^ contribution to thus contains both a constant part and an oscillatory 
part with frequency 2A?7. The spectrum is then given by the product of V^^'^ in 
Eq. (3.43) and in Eq. (3.51). 

Let us finally expand H, ei and €2 around the pivot scale k~^ = a{Xp)H{Xp), 
taken to be the logarithmic mean of the range of observable scales. At first order in 
the expansion, Vf'^ then reads 



2(1 + C)eip + Cesp + (2eip + eap) In 



(3.52) 



where the e^p and Hp are evaluated at kp. Assuming for simplicity that 0i and 02 
are ^-independent, the expansion of V^^^ around kp reads 

pose ^ ^2_^ 1^ _^2| |l-^2|l/2^ ^3_53) 

[cos {2k Ar] + 0i - 02) + n (2eip + eap) sin {2k Ar] + 0i - 6^2)] . 

At first order in slow roll, and in the limit <^ — 1, the modes v~ and are in 
the Bunch-Davies vacuum, and any modifications to induced by the bouncing 
phase disappear. This conclusion is also true for ? — )■ but this limit corresponds 
to a set of negative frequency modes v~ and (see Sec. 3.4). For <^ very different 
from or 1, oscillations can appear at leading order in the spectrum, rather than 
as small corrections. Note, finally that the form of obtained is reminiscent of 
the primordial spectrum obtained in trans-Planckian inflationary models but with 
some important differences. In trans-Planckian models [6, 7], the oscillatory term 
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has a In(fc) dependence rather than a dependence hnear in k. In addition, the 
amphtude of oscillations in trans-Planckian models is suppressed by the product 
of the slow roll parameters with the ratio i^mf/Mc ^ 1 where iJinf is the Hubble 
parameter during inflation and Mc is a new mass scale of order nip^. Here instead, 
the amplitude depends only on the deviation in amplitude from the vacuum state 
and is not suppressed by factors of order or by some new mass scale. 

4. Tensor perturbations 
4.1 Tensor perturbation variable 

In this section and the next, we compute the transfer of tensor perturbations through 
the bouncing phase and the primordial spectrum of gravitational waves. The evo- 
lution equation for the tensor modes h (where here h is not to be confused with 
the reduced Hubble parameter /i ~ 0.7 appearing in quantities corresponding to 
distances) reads 



where the subscript k on the mode functions h is implicit. Because tensor perturba- 
tions are not subject to gauge issues, h is everywhere well-defined and its evolution 
equation is regular. As for Vu{ri) in Sec. 3.2, we neglect the spatial curvature term 
and use the approximate form 



h" + [k^ - Vhiri)] h = with Vhir]) = + - 2/C 



(4.1) 



Vhiv)-^^ (2-ei) . 



(4.2) 



Using Eq. (3.15) the equation for h becomes 




(4.3) 



The solutions of this this equation before and after the bounce read 



h~ = .Jk^_ T{{k) Hl^J {kx.) + T:^{k) hI^J {kx.) , 



(4.4) 



= ^ T+ik) Hi'^ ikx+) + T+{k) Hi'^ ikx+) , 



(4.5) 



where, keeping terms up to first order in slow roll. 



3 



(4.6) 




and where the coefficients and are unknown functions of the wavenumber k. 
Matching the coefficients of /i"*" and h~ aX rj = 0, expanding the Hankel functions for 



\kx±\ <^ 1, and retaining terms in the expansion suppressed by negative powers of 
kArj up to linear order, one obtains 

T+ = (Tfe + i) e-*('=^'?-^^T), (4.7) 
T+ = Tf {Tk - i) e^(^^^— ^t), (4.8) 

where 

Contrary to the expressions for Ui and U2 in Eqs. (3.22) and (3.23) which were 
found to be exact at first order in slow roll, the truncated expressions for and 
T2' in Eqs. (4.7) and (4.8) are not exact at first order in slow roll since, in principle, 
higher powers t^, with p G N, contribute constant terms and terms of order ei to the 
expansion. At first order in slow roll, these terms read = (1 + 3pei/2)/[kAriy. 
These higher order terms {p > 1) can be safely neglected only if the constraint 
^ 0{ei) is satisfied. It can be verified that for Arj > 71 {At] = vr corresponding 
to a de Sitter bounce) and for any ei < 1 this constraint holds true for any n > 20 
[recall that k"^ = n{n + 2) for /C > 0]. 

Consistency with quantization again requires that the Wronskian condition be 
imposed on the tensor modes. As for the scalar modes, the Wronskian condition 
imposed prior and after the bounce reads 

\Tt\'-\Tt\' = T\k-\ (4.10) 
As for f , we now assume that T{ and T2 are parameterized as 



T{ = ^^1 A;-"T/2e^V'i and = —^2 k-^^'^e'"^^ . (4.11) 

Eq. (4.10) then imposes = = 1 and '01 = — 1- Again we drop the subscript 
on di so that di — > which gives 

T- = ^ A;-i/2e^^i ^nd T^" = ^ |l - k-^'^e'"^^ . (4.12) 



From Eqs. (4.7) and (4.8), the Wronskian condition on reads 

TT 

4 



\Tt\'-\n\' = -"rk^'{l + Tl). (4.13) 



Unlike for scalar modes, the Wronskian condition for /i+ only holds if <^ ^{^i) 

4.2 Primordial tensor perturbation spectrum 

For tensor modes, the primordial power spectrum "P^ reads 



k^ 647r 

V = 

^ 27r2 m2 
PI 



h+ 



T 

2 



(4.14) 
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Using the late-time expansion of the solution in Eq. (4.5), this becomes 



K = — 2 ikx^)''""^ 1^1^ - T,"-]' . (4.15) 

^ 27r2 a2 [r (1 - z/J sin (vrz/J]' ^ ii 21 v ; 

Once again splitting Eq. (4.15) into a standard part and a bounce-induced oscillatory 
part, = P'*'^ X we have 



16 1 

4 



pstd ^ ^ ^ ^ (kx+V-^"-^ , (4.16) 

- 2m>2[r(l-z.,)sin(7rz.,)f ^ 



= - A;|T+ -T2+I". (4.17) 



TT 



Expanding at first order in ei, one recovers, from 'P^^'^^ the standard result 



16 H'^ 



Vf = [1 - 2 (C + 1) ei - 2ei In {kx+)] , (4.18) 



vr 

PI 



while is analogous to the form of V'^^'^ with the substitution — > {i = 1, 2). 
Using Eqs. (4.12) and (4.6) and proceeding as in the scalar case, we find, up to first 
order in slow-roll, 

pose _ |^|2 ^ |l _ ^2| _ 2 |1 _ ^2|l/2 

X [cos(2/cA?7 + ^2) - Svrei sin(2A:A?7 ^1 - ^2)] • (4.19) 

As for the scalar perturbations, Eqs. (4.18) and (4.19) can then easily be expanded 
around the pivot scale k^. 

4.3 Tensor-to-scalar ratio 

The tensor-to-scalar-ratio r also splits into a standard and a non-standard part, 
^ _ ^std ^ j,osc_ rpj^g ratio of the standard parts of the tensor and scalar spectra yields 
the usual result at leading order, r^^'^ = p^td^-pstd _ iq^^^ while the non-standard 
parts of the spectra yield the additional ratio 



-pose \rp+ rp+\2 

' T _ -^2 
pose - \V+ -V2 



~ T>osc ~ |T/+ T/+I2' (4.20) 



Given the choice in Eqs. (3.34) and (4.11), this ratio depends on the values of q 
and on the phases ^1, 62, ipi and ip2, and it is scale-dependent through the 2k At] 
term in the sine and cosine arguments of V<^ and P^. The tensor-to-scalar ratio (4.20) 
therefore does not satisfy the standard consistency relation of single field inflationary 
models, for which one has r = — 8nT with nx = — 2ei the tensor spectral index. 
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5. The characteristic time scale Aij 



In previous sections we saw that the term 2kAri plays a crucial role in the modified 
spectra after the bounce, appearing in the arguments of the sine and cosine functions 
in both V: 



and V"'^. 

T 



We shall therefore now discuss the significance of the time 
scale Arj defined in Eq. (3.25). We first review the relevant field space dynamics as a 
function of the number of e-folds in slow roll infiation and then turn to the discussion 
on Ar]. 



5.1 Inflationary field space dynamics 

In infiation, the evolution of <y9 as a function of the number of e-folds is a known 
function, depending only on the parameters of the potential and on the initial con- 
dition for If at the onset of infiation. Indeed, as long as the slow-roll conditions are 
satisfied, i.e. when <^ 1, the horizon fiow functions introduced in Sec. 2.3 can be 
written in terms of field derivatives of the potential (2.23), 



ei 



1 fv_ 

2k 



V 



62 



K 



V 



V 



(5.1) 



which using the form of V{ip) in Eq. (2.23) gives 



ei 



2 ' 



£2 



8 (/i^ + 

K (/i^ — (/?2)^ 



(5.2) 



Thus, if the primordial spectrum is to be evaluated at some pivot scale /cp, the values 
of ei and €2 (and therefore that of the scalar spectral index rig) depend only on /i and 
on the value of the field, yPp, when modes of wavenumber fcp exit the Hubble radius. 
The end of the infiationary expansion is defined by ei = 1 at (/? = (yjo- This yields 



1+2^ 



1/2 



- 1 



(5.3) 



The slow-roll trajectory is then found from 



dip 



/i 



2 In 



(5.4) 



where cpi is the value of ip at some initial time. Inverting this expression, one can 
obtain the field value at Hubble exit of modes corresponding to the largest observable 
scales. 



-Wn 



fc 



exp 



8N, 
Kjj? 



(5.5) 
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where PVo(a;) is the principal branch of the Lambert function and ip^ is the field 
value evaluated at A^^, the number of e-folds between the Hubble exit of the largest 
observable modes and the end of inflation. Given a value of and a value of N^, 
(typically 40 < A''* < 60) ipi, is obtained from Eq. (5.5). This can then be used 
to determine the value of M by adjusting the amplitude of the spectrum of scalar 
perturbations to the COBE normalization. For simplicity, we neglect the minor 
modifications to the expression for the COBE normalization that are introduced by 
the new form of Vc_ and by fi/c ^ and simply write the normalization in the usual 
way, see e.g. ref. [58], 

-2 

-12 _4 



= 487ei 



xlO~'^m*. (5.6) 



Let us take ip-^ to be the field value at the onset of inflation (in standard inflation, 
9?i is given as an initial condition). The total number of e-folds of inflation A'lnf 
is obtained when ip^ of Eq. (5.3) is inserted into Eq. (5.4). It is useful, as a first 
approximation, to consider the case ip <^ fi for which A^inf can be approximated by 

A^inf = Ar((pe)^?/^'liif-l • (5-7) 



4 

Inflation should last at least 50 to 60 e-folds in order to solve the flatness, isotropy 
and horizon problems of Big Bang cosmology. Together with the value of ipe in 
Eq. (5.3), one may then use Eq. (5.7) to determine an upper bound for ip[ (in a small 
field model, the value of ip during inflation increases), which yields 

P>i < P>c exp . (5.8) 

This restricts ip\ to very small values. Recall once again that in standard inflation, 
ipi is given as an initial condition. We see from Eqs. (5.5) and (5.7) that the field 
evolution between (p^, and tp^ is independent of tpi but the total number of e-foldings 
of inflation is a function of ip\. 

5.2 Computing Ar^ 

Let us now turn to the case of the bouncing cosmology and the calculation of 
The conformal time interval Arj is obtained by integrating 

At] = ?7+ — ?7„ = / drj. (5.9) 

Jr]- 

Because the slow-roll contracting and expanding phases (sr) and the non-singular 
phase close to the bounce (b) cannot be approximated by a unique expression, it 
follows that for a symmetric bounce 

Af] = Af]- + Ar/b + Ar]+ = Ar^ + 2Ar],, = dr/b + 2 / dr],, . (5.10) 
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Here, ri~ and t]^ are the values of rj at which the cosmology respectively exits and 
enters the slow-roll regime in the contracting and expanding phases. The slow-roll 
parts of the integral (5.9), evaluated in the interval from r]^ to 77+ is given by 



dA^ 



dH{N) 



(5.11) 



where A^i = ln(ai/ab) and = ln(ac/ab). In these expressions, ai is the scale factor 
at the onset of inflation and is the scale factor at the end of inflation. Expanding 
In'H(A^) around A^p where A^p = ln(ap/ab) is the number of e-folds from the bounce 
up to the pivot scale, with Ni < Np < Nc, gives 



In-H(A^) = In-H(A^p) 



dlnH 



dA^ 



(iV-iVp) + 



1 d'lnH 



Np 



2 dA^2 



(N - Npf + 



Np 



-Ip 



){N-Np) 



^eipe2p(iV-iVp)' 



At linear order in e,-, one thus has 



niN) 



so that 



Since - A^n 



,{l-eip)Np 



1) 



'Hp(eip 

A''^, with 40 < A^^ < 60, we finally have 



,{l-eip){Np-Ni) 



Hp [1 



(5.12) 

(5.13) 
(5.14) 

(5.15) 



while from Eq. (5.4) one also has 
Np-N, = -fi' 



2 In 



(5.16) 



Here, the crucial point is that the field value at the onset of inflation ip[, although 
it is constrained through Eq. (5.8), is this time not given as an initial condition 
but depends on the detail of the bounce dynamics that preceed inflation and is 
an unknown function of at and /i. This is directly related to the fact that the 
exact dynamics from the bouncing phase intermediate between the deflationary and 
inflationary phases is not a known (or easily approximated) function of ab and /i. As 
a result, Ar/sr but also A?7b and the initial field value ipi, and finally Arj are quantities 
that cannot be obtained analytically. 

The value of At] can however easily be obtained by integrating the background 
cosmology numerically for any value of ab and /i. We do so in Fig. 1 where we plot, 
in the full lines, the quantity aoA?7, which is the theoretical frequency of oscillations 
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for physical wave numbers, for the renormahzed potential parameter /i in the range 
2.0 < /i < 4.5, over the allowed range for at [see Eq. (2.24)] and in the range 
10~^ < \n,c\ < 10^ 

Furthermore, remembering that in the case of a /C = +1 universe, the wavenum- 
ber k is discrete with k = (n + 2) for integer n, the arguments of the sine and 
cosine functions in Eq. (3.53), which are relevant for the oscillatory features in the 
primordial spectra, read 

2A;Ar/ ~ 2 1^72 + 1 - (mvr + 5r/) (n > 1). (5.17) 

Here, m is an integer and hence < 6r] < it, such that, in Eq. (3.53), one can make 
the substitution 

+ 1\ 

2k At] 2 ao67] = kphysi2ao67]) (5.18) 

V '^0 / 

for large n, where fcphys, as previously, is the physical wavenumber and is the scale 
factor today. The effective frequency is therefore 2ao6ri rather than 2aoA?7. The 
quantity aoSr] is shown in the circles of Fig. 1. The quantity Srj itself is shown in 
a two-dimensional plot in Fig. 2 in the ranges 3.0 < /i < 4.0 and 2.04 < T < 2.08 
and for 10~^ < \^k,\ < 10'^. This reduced range for fi corresponds to Ug — 0.95, in 
agreement with the observationally constrained spectral tilt of the primordial power 
spectrum of scalar perturbations, while the range for T was chosen in order to satisfy 
observational bounds on Q/^. 

Fig. 1 demonstrates that for fi = 3.5 and /i = 3.0 the effective frequency drops 
significantly in some narrow bands of T [or equivalently bands of Ob, see Eq. (2.24)] 
while Fig. 2 shows the existence of this band for the case 3.0 < ^ < 4.0 in more 
detail. These bands correspond to values of Arj ~ 2tt with drj <^1. The existence of 
these narrow bands is crucial. This is because, given that oq is at the very least equal 
to 20/i~^ Gpc, the effective frequency of oscillations, ao^ry, can only be small enough 
for such features to be observable in the CMB or in the matter power spectrum if 
5rj < 1. 

In fact, the results of Figs. 1 and 2 suggest a new possibility for an indirect deter- 
mination of the spatial curvature fi/c should oscillations linear in k and attributable 
to a bouncing cosmology of the type described in this paper be measured in the 
temperature anisotropy map and/or in the matter density power spectrum. Given 
accurate measurements of and ao(5?7, a good estimate of can be obtained through 
Eqs. (5.1) and (5.2). Given an accurate measurement of ao5?7, the value of T is known 
so that the value of the scale factor at the bounce is recovered. Furthermore, given 
V{ip), the value of /i, say of order 3.0, and the value of Cb are enough to obtain A?7 
numerically. Finally, given the knowledge of Arj and that of ao5?7, the value of Qq is 
easily obtained. Combined with Hq, it yields fi^:- 
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At last, one can estimate the number of oscillations per decade D in the wavenum- 
ber k from 

Aro3e(D) = ^10^(^), (5.19) 

where i{D) is an integer. For example, for yU = 3 nip^ [see Eq. (2.23)], Cb — 2.5 x lO^fpj 
and setting 6r] = 0.001, the observable range 10~^ hMpc~^ < /cphys < 10/iMpc~^ is 
10^ < < 10^ and \Qjc\ < 10"^ This gives a total of four decades, D = 1,2,3,4, 
with i{D) = 3,4,5,6. In this case, A/'osc (2)) ranges from 0.3 (for D = 1) to 3 x 10^ 
(for D = 4). If instead, 5r] = 0.01, then A/'osc (2)) ranges from 3 (for D = 1) to 3 x 10^ 
(for 2) = 4). To see that these estimates are correct, see Sec. 6.2. 



\i = 4 n = 3.5 n = 3 n = 2.5 \i = 2 




Y Y Y Y Y 



Figure 1: Numerical integration of agAr/ and ao5r]. The full lines are aoArj, while the red 
circles are aodrj. For /i > 3, "/ig > 0.95 and Arj ~ 27r so that there exist regions in which 
a()6r] drops to significantly lower values, see Fig. 2. 



6. CMB and matter power spectra 

In previous sections, we derived the scalar and tensor primordial power spectra and 
showed that oscillations induced by the combined effect of the deviation of incoming 
perturbations from the vacuum state and by the characteristic time scale Arj appear 
in and Vh- In addition, we found that measuring the scalar spectral index Ug 
and the frequency of oscillations makes possible, within the scope of this model, an 
indirect determination of the values of and of Q/c today. 

The oscillations can however only be measured if they are transferred to the 
angular power spectrum of the CMB (the C^'s) and/or to the matter power spectrum 
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Figure 2: Numerical integration of Sij. The thin lines represent isohnes on a loga- 
rithmic scale, the red one denoting \0,ic\ = 2 x 10~^. In the range 3 < ;U < 4, ris ~ 0.95. 

Vs- The transfer function relating Vs to V( is simply a fc-dependent multiplicative 
function and oscillations can be expected to be transfered from to Vs unchanged, 
see Sec. 6.2. The Ces, on the other hand, are related to V,^ through a complicated 
expression so that the transfer of oscillations to the angular power spectrum of the 
CMB is by no means obvious. This question is investigated in this section. We 
first write down the analytic expressions for the angular power spectra in the range 
1 < £ < 10. We then give numerical examples of Q over the entire i range using the 
CAMB code [54]. 

6.1 Scalar anisotropy spectrum 

For \Qk:\ ^ 1 and for £ < 10 but for K, = +1, the scalar multipole moments C| can be 
related to the primordial scalar curvature perturbation spectrum V( by [68, 69, 70] 




2 




n=2 
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with 

f ^ 

- r? for £ < n 




for £ > n 



and where P_\^il^jX associated Legendre polynomials while X = r\Q — r\\^^ ~ 

^h/o-q = \/Wk\ is the (conformal) radial distance to the last scattering surface, with 
rjo and r^igg, the conformal times today and at the time of last scattering, respectively. 
With the exception that the integral over k should be replaced by a sum over n, the 
textbook expression used in the case /C = 0, see [70, 71], is a very good approximation 
to Eq. (6.1) for \Qic\ <^ 1. The expression for the scalar C^'s can be carried over to 
the case of the C/s for tensor perturbations. 

The presence of the slowly oscillating sinusoidal functions in introduces an 
oscillation in the kernel of the sum in Eq. (6.1). In particular, one can in general ex- 
pect constructive or destructive interference between the oscillations of the -P_i/2+n'^ 
having as a characteristic scale x with the oscillations of the cosine and sine terms 
in V(; of characteristic frequency 2Ari. Summed over n, this induces slow oscillations 
in the C^'s at small i values. The precise shape and locations of the deviations from 
the standard C^'s induced by the oscillations depend on fljc via x ci-nd on Srj while 
their amplitude depends on q. 

The COBE normalization can be expected to be modified in two ways. First, 
note that for ^ > 1 the constant term in V^^^ is equal to 2^^ — 1. The amplitude 
of the spectrum is therefore modified by this first factor. Secondly, the amplitude 
of the spectrum will also be modified by a contribution coming from the oscillatory 
part. This contribution depends, once again, on on x and on Ar]. 

In Figure 3, we present the results of a numerical simulation using the CAMB 
code [54]. The quantities shown in the dashed red and full blue lines are defined by 

AC, = C™^o^^-C™ (6.2) 

where the theoretical Q's are the C^s obtained, respectively, either from an inflation- 
ary cosmology or from the bouncing cosmology described in this paper. The error 
bars are the WMAP error bars for the measured C^'s. The AC^'s were computed for 
^ = 1.01, 1.1 and 1.2 and for 6t] = 0.001 and 0.01. These results were obtained for 
/i = 3, iV^ = 50, T ~ 2.055. For this choice of parameters, ei = 0.0025, eg = 0.045, 
A'^inf = 65, and for a period of radiation and matter domination (neglecting reheating) 
lasting 65 e-folds [72], such that Q/c = —0.002. The figure mainly illustrates that at 
the frequency 6ri = 0.01, and for q = 1.01, the oscillatory features in the primordial 
spectrum appear to improve the fit to the WMAP data in comparison with what 
is obtained in the absence of oscillations in P^. For larger values of ^, such as 1.1 
and 1.2, the effect of the oscillations worsens the fit near the first accoustic peak at 
i ^ 200 over a widening range of i values as <^ is increased. At higher frequency, for 



-28- 



6ri = 0.01, the resulting Q's have large fluctuations that appear to be phase-shifted 
with respect to the result obtained in the inflationary case. At low i, oscillations 
are not visible but there remains a change in amplitude with respect to the AC^'s 
obtained from inflation. This could be expected from Eq. (6.1) for the values of q 
and St] chosen. 




10 20 50 100 200 500 1000 10 20 50 100 200 500 1000 



Figure 3: ACi spectra for ? = 1.01, 1.1 and 1.2, from top to bottom respectively and for 
6r] = 0.001 and 6r] = 0.01 from left to right, generated using CAMB. The error bars are the 
WMAP error bars on the measured values of the C^'s, the red dashed lines are the AC^'s 
for an inflationary cosmology and the blue line in full represents the AQ's for the bouncing 
cosmology presented in this work. 

6.2 Matter density power spectrum 

The matter power spectrum Vs is related to the primordial power spectrum P({k) 
by a /c-dependent transfer function T{k), i.e. one can write Ps{k) = T{k) Pc^{k). 
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Roughly speaking, T{k) ~ 1 for small k, and T{k) oc at large wavenumbers. The 
oscillating behaviour observed in the primordial perturbation spectrum is therefore 
expected to be transmitted to the matter power spectrum. Again, using the CAMB 
code [54] with standard cosmological parameter values, we have computed the theo- 
retical matter power spectra for the values of the parameters specified in Section 6.1. 
The results are shown in blue in Fig. 4. The theoretical spectra were then convolved 
with the observational window functions of the Sloan Digital Sky Survey (SDSS) [73] 
and are shown by the red dots in Fig. 4. These theoretical results are compared with 
the SDSS data, shown in black in the figure. The convolution evidently smoothes 
out the oscillations, and it is clear from Fig. 4 that there exist values of <; and Ar] 
such that the resulting matter spectra are fully degenerate with those obtained from 
a standard slow roll power spectrum. The range of scales displayed is centered on 
the second decade {D = 2) of observable scales. The number of oscillations in this 
decade is about 3 or 4 for 61] = 0.001 and about 20 to 30 for 61] = 0.01. These are 
the estimates that were found in Sec. 5. 

7. Conclusions 

This paper refines and completes the analysis begun in Ref. [38] and continued in 
Ref . [36] . In Ref. [38] , the authors focused on the immediate vicinity of the bouncing 
phase and conducted a detailed analysis of the bounce-inducing background cosmol- 
ogy and on the transfer of fluctuations through the bounce by modifying the exact 
solution obtained when considering a de Sitter universe with closed spatial sections. 
In Ref. [36] , the framework used in Ref. [38] was used to show that for symmetric or 
quasi-symmetric general relativistic bouncing cosmologies with /C = +1 the peak in 
the potential for the variable u could never be large so that the metric perturbation is 
always sub-Hubble in the vicinity of the bounce. A cosmology smoothly connecting a 
contracting phase and a bouncing phase to an inflationary phase was then proposed 
and analyzed at both the background and linear perturbations level. In the present 
work, we have provided a detailed summary of the background cosmology exploited 
in Ref. [36] and have performed a much more detailed calculation of the transfer 
of perturbations through the contracting, bouncing and inflationary phases of the 
cosmological background discussed in Ref. [36]. We have parameterized the initial 
state of perturbations prior to the bounce in terms of the Bunch Davies vacuum state 
and used the solution for the Mukhanov-Sasaki variable (far from the high spatial 
curvature region) in terms of Hankel functions, to reduce the number of unknown 
parameters needed to define the initial state of first order perturbations, see Sec. 3.4. 
We then computed the effects of the bounce and the choice of initial conditions on the 
scalar and tensor mode primordial spectra "P^ and Vh- We discussed the modified low 
multipoles of the CMB angular power spectrum, the modified COBE normalization 
and also provided numerical evidence, using CAMB, that the C^'s and matter power 
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Figure 4: Matter power spectrum for = 1.01, 1.1 and 1.2, from top to bottom respec- 
tively and for 6r] = 0.001 and 5rj = 0.01 from left to right. The black dots with error bars 
are the data points from the Sloan Digital Sky Survey (SDSS) [73], the blue line is the 
theoretical prediction of the bouncing cosmological model obtained using CAMB, and the 
red dots are the simulated data points obtained by convolving the blue line with the SDSS 
window functions. 



spectrum Vs are affected by the combined effect of the choice of initial conditions 
and the bouncing cosmology. 

The crucial point is the appearance of oscillations in the power spectra. These 
oscillations mainly depend on the initial state of the perturbations through the free 
parameter q and on a new cosmological scale Arj. The former sets the amplitude of 
the oscillations while the latter sets the oscillatory frequency. As shown in Sec. 6, 
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there exists values of ^ and Ai] for which the oscillations induced by the bouncing 
phase and by the choice of initial conditions do not conflict with the WMAP and 
SDSS data. In fact, there appears to be parameter ranges in which the C^'s and 
the Vs derived from V( are fully degenerate with those obtained from a standard 
slow-roll primordial spectrum. The tensor-to-scalar ratio is also modified. It departs 
from the standard result in two ways, first by a modification of the overall amplitude 
and secondly by the existence of a scale-depedent oscillation, see Eq. (4.20). 

We also identified a new way of indirectly measuring the spatial curvature of 
the universe, in Sec. 5, assuming the oscillations can be attributed to a bouncing 
cosmology of the kind described in this paper and can indeed be measured. The 
scale factor at the bounce as well as the model parameter n determine both Q/ck and 
A?7, thereby establishing a one-to-one relationship between the two. If the frequency 
of oscillations could be measured together with the spectral index of P^, then Q/c is 
determined, as well as a^- 

In Appendix A, we discussed the issue of the growth of perturbations in the 
contracting phase and pointed out the possible necessity of studying asymmetric 
bouncing cosmologies. This will be the subject of future work. Finally, a well-known 
problem of contracting spacetimes with spatial curvature is the instability of the 
background evolution [36, 37, 52]. Again, this will be the subject of future work [74]. 
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A. Time evolution of curvature perturbations 

In this Appendix, we discuss the time evolution of curvature perturbations. We 
first provide some numerical examples for symmetric as well as for non-symmetric 
bounces and then discuss the growth of perturbations in the contracting phase. 

A.l Numerical examples 

Before gauge fixing, the scalar part of the perturbed metric is given by [59] 



ds^ = {r]) [- (1 + 2A) dr]^ + (1 + 2C) -f^jdx'dx^] , (A.l) 



where A and C are small but gauge-dependent scalar perturbation variables. The 
gauge-invariant scalar degrees of freedom are [59] 



^ = A + H{B -E') + {B - E')' , 
m = -C -U{B - E') 



(A.2) 
(A.3) 



where B and E are the scalar parts of the g^i and g-ij components of the perturbed 
metric and where $ and ^ are equal in the absence of anisotropic stress in the stress- 
energy tensor. The curvature perturbation on comoving and on constant energy 
density hypersurfaces, denoted by TZ and C, respectively, are given by [75, 76] 



n = c-H 



p + p 



and 



C = c-n 



5p 
P' 



(A.4) 



where q is the scalar potential of the momentum density qi of the perturbed energy- 
momentum tensor, defined as qi = Viq + with V^g* = 0. Both TZ and ( can be 
expressed in terms of gauge-invariant quantities as 



7e = ^ + 



2n ^' + 
3{h^ + k:) 1 + w 



(A.5) 



and 



C = <i'- 



2n 



3{1 + w){H^ + IC) 



1 / k 
3 



(A.6) 



On super-Hubble scales and in the absence of isocurvature perturbations, C' = if 
$' + "H^ does not diverge (this can be seen from the equation of motion for $), and 
TZ ~ — C provided /C ~ 0. Hence, both ( and TZ are conserved after Hubble horizon 
crossing. In the absence of isocurvature perturbations, this holds in very general 
circumstances. These two quantities are therefore of primary cosmological interest. 
For scalar field matter, TZ can also be written in Newtonian gauge as 



7^ 



6ip 

V 



(A.7) 



and it can be seen from (3.8), taking /C = 0, that TZ is related to v hj v = —zTZ. 
If /C/a^ ^ k"^, we may work with v and z rather than with v and z. In this case, 
C ~ —TZ = v/z, and knowledge of v is directly usable to obtain (. 

In Figure 5, we show the results of a numerical integration of the background 
cosmology, see (2.11) to (2.12), and curvature perturbation on equal density hyper- 
surfaces (, by integrating (3.4) and then using (3.3) and (A.6), for three sets of 
initial conditions at the bounce with which the bouncing phase gets progressively 
more asymmetric. 
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A. 2 Growth of perturbations in the contracting phase 

We now discuss the issue of the growth of perturbations in contracting spacetimes as 
discussed in Ref. [65], see also [67]. Let us consider the simple example of power- law 
contraction of the universe in which a{r]) ~ where p = 2/3(1 + w) is positive 

and X = Irj — T]c\ for some characteristic conformal time rjc < 0- li p < 1, t] < r]c and 
if p > 1, ?7 > ?7c so that for any p > the universe contracts as rj goes from — oo to 0. 
In the power law case, it is easy to solve the evolution equations for either $ or TZ. 
The solutions can be expressed, on super-Hubble scales, as the sum of two modes, 
one constant and one time-dependent, 

$fe ~ + C$,2 (A.8) 

and 

TZk ~ C7^,lX(3p-^)/(^'-^) + Cn,2 • (A.9) 

When ?7 — )■ as the bounce is approched, in a matter- or radiation-dominated con- 
tracting phase, p = 2/3 oi p = 1/2, the C$,i mode diverges as and respec- 
tively, while the C-ji^i mode diverges as X"^ and X~^ respectively. If p ~ as in the 
ekpyrotic scenario [32, 33], the C$,i mode grows as X"^"^^ but the Cti,i decays as 
^i-2p_ If y; < —1/3, corresponding to a > 0, then p > 1 and X grows. For w ~ — 1, 
p ^ 1 and the C$,i and C-ji^i modes grow as and respectively. How- 

ever, in this case, X ~ 1/a and a ~ e^*, so that the conformal time interval that a 
mode remains super-Hubble in this regime is exponentially small so that both modes 
remain approximately constant. 

It is clear that the perturbation variables $fc and TZk do not evolve with the same 
power in X. Furthermore, they are related to each other through their derivatives, 
so that, generally, the time-dependent mode of $ does not correspond to the time- 
dependent mode TZ, and vice versa. The example given above also implies that in 
a given cosmological setup, perturbation theory may be preserved or the possibility 
that it be violated be alleviated provided the right choice of perturbation variable 
is made [65, 66]. In an explicit model where modes evolve in and out of the Hubble 
radius, and where the equation of state parameter varies, as is the case for scalar field 
matter, the issue of whether modes grow large or not evidently becomes more subtle. 
Let us ask which variables must remain small. In longitudinal gauge, the scalar part 
of the perturbed metric is diagonal, B = E = 0, so that the line element takes the 
form given in (3.1). It is legitimate to ask whether and $ should remain small, or 
whether it is A and C that should remain small. If the latter choice is made, then the 
gauge must be suitably chosen so that A and C do remain small in the cosmological 
evolution being studied. As in the example above, the longitudinal gauge may not be 
appropriate if $ is a growing function of conformal time. In order to decide whether 
there exists a gauge in which perturbation theory remains valid, it is interesting to 
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look at the intrinsic curvature of spatial sections and at the extrinsic curvature of 
comoving surfaces. The intrisic curvature of spatial sections in the comoving gauge 
is given by 

(3)i? = ^ + l(A + 3/C)7^, (A.IO) 

while the scalar part of the extrinsic curvature on comoving hypersurfaces is given 
by 



K' = - (1 + 27^) 5] + 



an ^ 2 A 



(A.ll) 



where A~^ is the inverse Laplacian. The evolution of each term in ^^^R and Kj 
as a function of X is shown in Table 1. In all cases but the ekpyrotic scenario, 
perturbative terms to ^^'^R and Kj grow faster than background terms. In order 
for perturbation theory to be preserved, the contracting phase must therefore be 
short - this points to the study of strongly asymmetric bouncing cosmologies [67] 
- or the contracting phase must be induced by a potential for (p giving rise to very 
large w; see however [77]. In this work, see Sec. 3.4, we assume that the deflationary 
phase is short enough that perturbations remain small in the rapid contracting phase. 
Figure 5 shows that this is in fact a realistic assumption. We further assume that 
the amplitude of fluctuations at the onset of the deflationary phase is small. We 
thus choose initial conditions for the Mukhanov- Sasaki variable v with an amplitude 
deviating only slightly from that of vacuum fluctuations, see Sec. 3.4. 



Type 


P 


X 




n/a 




nn/a 


TZ/aU 


TZ' 


Ekpyrotic 


< 1 


\ 








x-^p 


X^-sp 


x-^p 


Radiation 


1/2 


\ 














Matter 


2/3 


\ 










X-4 




Exponential 


> 1 


/ 


X2+2/P 






X3+3/P 


X^+s/p 


X2+2/P 



Table 1: Super-Hubble conformal time evolution for each term in the intrinsic and extrinsic 
curvatures in the comoving gauge. 
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Figure 5: From left to right, number of e-foldings A^, Hubble parameter H, potential 
Vu [true potential (blue), approximate potential (red dashed)] in rows 1, 3 and 5. The 
normalized curvature perturbation C/Ci [real part (blue), imaginary part (red dashed)] for 
n = 5, 10, 15 in rows 2, 4 and 6, with Q = C,{ti) the curvature perturbation at time tj. 
From top to bottom : T = 2.05, = (top rows), T = 2.05, = — ?Tipj/6 (center rows), 
and T = 2.5, (^b = —rri-p^/Q (bottom rows). 
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